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1. Introduction 

Over the decades, polylogarithms have gained importance in perturbative Quantum Field The- 
ory ever since the first occurrences of the dilogarithm in early QED (e.g. [|30|]) and in results for 
one-loop integrals. As the complexity of the computations has increased, the literature on Feyn- 
man integrals has gradually absorbed a variety of generalizations, including Nielsen and classical 
polylogarithms [ |26| , 28, ^3|], harmonic polylogarithms [ j29| ] and later generalizations Jl8| , f3[ |], g], 



(some of which were previously known to mathematicians by the name of hyperlogarithms p5[]), 



and variations on multiple polylogarithms [n, 19]. For several approaches to computing Feynman 



integrals, it is useful to represent these functions as iterated integrals. 

In this talk we discuss a class of functions which was studied in reference [^]. These functions 
are closely related to (and contain) the multiple polylogarithms of Goncharov [ |l9"l ] and admit spe- 
cial properties which are useful for the computation of Feynman integrals. In section ^| we briefly 



review an important result by Chen Q14] ] on the conditions for an iterated integral to give a well- 
defined function of several variables. In section [| we define a map whose image satisfies these 
properties and use it for the construction of the mentioned class of functions. Section |] briefly 
shows how these functions can be used in a systematic approach for integrating over Feynman pa- 
rameters, presented in reference [Sh. Algorithms and a computer program for using this class of 
functions will be the content of another publication, which is currently in preparation. 



2. Homotopy invariance and Chen's theorem 

The property of homotopy invariance is best discussed when viewing iterated integrals as 
integrals along paths. Let k be the field of either the real or complex numbers and M a smooth 
manifold over k. Let a piecewise smooth path y on M be given by a map y : [0, 1] — > M. Two 
such paths /i, 72 are said to be homotopic if their endpoints coincide: 71 (0) = 72(0) = xq and 
7i(l) = 72(1) = X\, and if furthermore one path can be continuously transformed into the other. 

Let COi, (O n be smooth differential 1 -forms on M and let us write 7* (to,) = ft(t)dt for the 
pull-back of each 1-form to the interval [0, 1]. The iterated integral of COi, ft),, along 7is defined 
by 

(O n ...(O l = fn{tn)dt n ...fl(tl)dh. (2.1) 

Jy Jo<h<...<t„<i 

We will use the term iterated integral for ^-linear combinations of such integrals. 

To give an example, the multiple polylogarithms in one variable can be written as 

U, 1 ,...,, v (z) = (-l) r / (2-2) 
Jy 

where 7 is a smooth path in C\ {0, 1} with endpoint 7(1) = z and the sequence of 1-forms is built 
up from the set Q.[ = {(Oq, a>i} with (Oo = y, (Qi = j^. Note that if we extend the set Q[ only by 
1-forms of the form f dt, where / is a rational function, possibly involving further parameters, and 
use it to define a new class of iterated integrals in a similar way, then these functions will still be 
defined on a one-dimensional space, corresponding to the one endpoint variable z. In this case let 
us speak of iterated integrals in one variable. The manifold M is an open subset of C. In the next 
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section, M is a certain open subset of C'\ and we use 1-forms with several dt\ , dt n to construct 
a class of functions in n variables z\, ■■■,z n , given by the coordinates of the endpoint of a path. 



The iterated integral in eq. is a meaningful expression for multiple polylogarithms because 
it depends locally on the endpoint variable z but not on the path 7. In fact if in eq. we replace 7 
by another path homotopic to 7 then we obtain the same function. This property is called homotopy 
invariance. In the case of one-fold integrals of a 1-form ft) one can show that homotopy invariance, 

co = co for Yi , 72 homotopic, (2.3) 
■hi '~h 

is true, if and only if co is closed. 

For iterated integrals the condition is more complicated. It was studied in a very general setting 



in Chen's foundational work on iterated integrals [|14j] and we want to briefly rephrase the statement 



which is relevant in our context. To this end we consider tensor products of differential 1-forms 
ft)i <g> ... <g) co m over some field K (Ik (which will typically be the field of rationals Q in the sequel) 
for which we use the customary bar notation [c0\ \...\co m ]. 

Let £2 be a finite set of smooth 1-forms on M and let D be the ^-linear map from tensor 
products of such 1-forms to tensor products of all forms on M, defined by 

m m— 1 

D([co x \...\co m }) = ^[fti 1 |...|fti,-_ 1 |Jfti,-|fti,- +1 |...|ft) m ]+ £ [fti 1 |...|fti,-_ 1 |fti,-Afti,- +1 |...|ft) m ]. (2.4) 



1=1 



We furthermore define 



B m (Q) = < t, = ]T £ Ci u ...fy [ft);, |...|ft),,] with £ K and ft), £ D. such that Dt, = > (2.5) 

I l=0h,-,il J 

which is a vector space over K. We call the elements of this vector space integrable words (or 

bar elements) in £2 and the equation Dt, = is known as the integrability condition. Now on 

the elements of B m (Q.) let us consider the integration map, defined by simply integrating over the 



1-forms according to definition fl2.1p 



L L c h,..j, W-K] ^ L L c h,.-,h / °v -®h ( 2 - 6 ) 

2=0 11,..., 11 2=oii,...,ii J y 

Chen's theorem now states, under some conditions on Q. which will always be satisfied in the 
sequel, that this integration map gives an isomorphism from B m (Q.) to the set of homotopy invariant 
iterated integrals in the set of 1-forms in Q. of length less than or equal to m. 

In other words, when we apply the integration map to an integrable word, we obtain a ho- 
motopy invariant integral. The reverse is also true: any linear combination of tensor products, 
corresponding to a homotopy invariant iterated integral in the above sense, is an integrable word. 
In the following we implicitly use this isomorphism and represent a homotopy invariant iterated 
integral by its bar element (this requires fixing basepoints: we shall demand that the regularised 
value of all functions at the origin is zero). Let us consider all integrable words of a given alphabet 
£2 up to length m, B m {Q.) and define S$ m (CI) to be the ^-vector space of all the corresponding 
homotopy invariant iterated integrals, obtained from these words via the integration map of eq. 



2.6. For notational convenience, the bar notation, used for words in B m (Q.) above, will from here 
on denote the corresponding functions in ,% m (£2) as well. In the following section we explicitly 
construct this vector-space of functions for a specific choice of £l. 
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3. Universal multiple polylogarithms in several variables 

From now on, let K = Q. Extending the set Cli of the previous section, let us define the 
auxiliary set of differential 1 -forms 

- _ \ dt v dh hdh (n? =2 r,-)<fti 1 

n ~\h '/i-r*ift-r-' uuu-i j { } 

on an open subset of C with coordinate t\. As the 1 -forms in Q. n are closed and furthermore the 
wedge product of each pair of 1 -forms in Cl n is zero, the integrability condition is trivially satisfied 
for any word in the letters of Cl n , and therefore B m (£1„) is spanned by the set of all words of length 
< m in Cl n . Therefore it is trivial to obtain the homotopy invariant iterated integrals 38 m (fl„) . 

However, all of these are homotopy invariant functions of only one variable, the endpoint of 
the path corresponding to the integration over dt\. By a slight abuse of notation let us call this 
variable t\. The other parameters t 2 , t n in the 1 -forms have to be treated as constants up to now. 

Let us now consider the set 

( dty dt n d(U a<i<h ti) , a , 

&n = \ — , -, — , „ ~ ~ : where 1 < a < b < n } . (3.2) 

L h tn Ha<i<bU- 1 J 

For example, in the case of three variables we have 

_{dt\ dt 2 dt 3 t l dt 2 + t 2 dt i t 2 dt 3 +t 3 dt 2 t\t 2 dt 3 + t 2 t 3 dt\ +t]t 3 dt 2 } 
3 ~\t^ t 2 ' t 3 ' ht 2 -\ ' t 2 t 3 -l ' ht 2 t 3 -\ )' 

Note that for n > 1, in contrast to the previous case, not every possible word in Cl n belongs to 
B m (H„). We explicitly construct SS m (H„) by a map 

(«n) , 0.4) 

which is defined as follows. 

Let F n be the vector space of rational functions of t\, t n with coefficients in Q whose de- 
nominators are products of elements in the set {t\,.. . ,t n ,Y[ a <i<bti — 1}, for 1 < a < ft < «. In the 
following we write iterated integrals in SS m (£2„) as \g\dt\ \g 2 dt\ \...\g n dt^\, where gi € F n . Differ- 
entiation of the iterated integrals in ,% m (£l n ) with respect to t\ can be computed by 

d d 

-^-[gidh]=gi and — [gidh \g 2 dt\ \...\g n dh] = gy [g 2 dt\ \...\g n dt\] for n > 2. (3.5) 

The map \ff will not change the differential behaviour with respect to t\ , so we impose the differen- 
tial equations 

d d 

— y{[gidh])=g l and —\\f{[g l dh\g 2 dti\...\g n dh}) =g\V(\gidh\...\g n dh\) for«>2. 

Differentiation with respect to t 2 is not defined on § 6 g$ m [fl n \ but we want y(^) to have a 
well-defined differential behaviour with respect to t 2 . To this end, we consider an auxiliary operator 

(a n )^F 2 ®3§ m (a n ) (3.6) 
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where d tl is defined by the following properties: 

(a) On rational functions it acts as differentiation with respect to t 2 : d t2 g = -^-g. 

(b) It commutes with differentiation with respect to ti: d^j^E, = -^d t2 t,. 
Using property (b) we obtain 

d 



dti 



d, 2 [g i dt l \g 2 dti\...\g n dt l 



and therefore 



d, 2 [gidti\g 2 dt l \...\g n dt l 



d t2 -^[gidt 1 \g 2 dti\...\g n dt l } 
d t2 gi [g 2 dti\...\g n dt\\ 



dt[d t2 gi [g 2 dt[\...\g„dt[ | 



(3.7) 
(3.8) 

(3.9) 



Note that on the right hand side of the last equation, d t2 acts on an iterated integral of length n — l, 
so we have a recursive procedure to compute d, 2 E, , the last stage of the recursion given by property 
(a). In (|J!9| ), the integral is computed by decomposing d t2 g\ E F n into partial fractions with respect 
to t[, and using integration by parts and the formula ( |3.5[ ). 
Now we define y such that it satisfies 

d 



dt, 



and analogous differential equations with respect to the remaining parameters f,. The constants are 
fixed by demanding that y(^) has a finite expansion at the origin of the form 



I A,...,,„(log fl )"...(log ?n ) ! " 

0<h,...,i n <N 



(3.10) 



where fa , n is analytic at the origin and vanishes at the point t\ = ... = t n = 0. Together with the 
above differential equations this determines the map y. This map is closely related to constructions 
which were recently introduced to the physics literature as the 'symbol' p0| , p| , 16]. 
To give an example, we apply y to 



dti 
1^1 



t%dt\ 



hh-l 



fn, 



We obtain 



1 



dt 



i 

dt 2 



ti dt 2 + t 2 dt 



dti 
dti 



1 



+ 



ti-l 



+ 



1 

h h 
dt\ 



1 

dt 2 



ti dt 2 + 1 2 dt\ 



7 + 

?1-1 t 2 f 2 — 1 

This iterated integral y(^) is in fact equal to the multiple poly logarithm in two variables 
Lit, i (?i , h) while the expression t, coincides with this function only on a one-dimensional subspace 
for fixed t 2 and does not capture its differential behaviour with respect to t 2 . 

Let 38 (0„) = Y.m>o^m (n„) denote the vector space of integrable words of all lengths m > 0. 
It was extensively studied by one of us in the context of the moduli space of curves of genus zero 
with m + 3 marked points, ^#o,m+3 , in reference [J7|] and we refer to this work for details and proofs 
of the following properties: 



dt 2 



1 



ti dt 2 + 1 2 dti 



hh - 1 



(3.11) 

(3.12) 
(3.13) 
(3.14) 
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• The elements of 3S(Q. n ) are homotopy invariant and therefore they are functions of n vari- 
ables. 

• 38 (n„) contains the multiple polylogarithms of Goncharov. 

• There is an explicit basis for 38(Q. n ) in terms of the map *P. There is a decomposition 
38 (Q n ) = *¥(38{Q. n ))® ...® x V{3§(Q.2))®^(^{^-\)) and we have an explicit basis for 
each 38 (ClA given by the set of words in by the discussion after (fO|). 

• 38 is closed under taking primitives. 

• The limits of elements of 38 (£l„) at t„ equal to and 1 are «2°-linear combinations of elements 
of 38 (Q„_i), where 0£ is the Q- vector space of multiple zeta values. 

As a consequence of the latter properties, we can evaluate definite integrals of the type 

fdt^fjPj, pj£SS m {Q n ), fj€F n . (3.15) 
Jo j 

The result will be a iT-linear combination of elements of 38 m (Q. n -\) multiplied by elements of 
F n -\. This can be iterated. Our forthcoming publication will feature algorithms and a computer 
program for the computation of such integrals. 

4. Application to Feynman parametric integrals 



The integrals ( |3.15 ) play a role in pure mathematics, such as in reference [Q, [l0| ], and in 



physics in the context of deformation quantization [|l7|], superstring theory [31], Schnetz' model of 



graphical functions [ p2[ ] and in perturbative quantum field theory. Here we focus on the latter and 
give a very brief outlook on how the use of 38 m (O n ) can facilitate the computation of Feynman 
integrals. We follow the approach of reference []8|] which in combination with the use of hyper- 
logarithms already led to new results for certain integrals relevant in QCD [[]]]. Other results were 
recently obtained by similar strategies of integrating over Feynman parameters, e.g. in [Jl3||. 

To begin with, we consider a primitive (subdivergence-free) overall logarithmically divergent 
vacuum Feynman graph, giving rise to a finite integral 



llUdxi )S \l-Vxi\-spj (4.1) 




over N Feynman parameters with being the first Symanzik polynomial (see e.g. [|27|, [220). Ref- 
erence ^ provides a polynomial reduction algorithm which for any ordered sequence of the N 
Feynman parameters X = (x CTl , ...,x Gn ) gives a sequence S*x = (Si, Sn) of sets of polynomials 
in the Feynman parameters. Without repeating the details of this algorithm here, we recall that one 
can evaluate / by iteratively integrating over the Feynman parameters in the order A if in each set 
Si € 5P\ all polynomials are linear in the corresponding parameter x 0j . 

Note that the latter is exactly the condition for the integral to be computable by our program, 
using the functions in 38 m . Indeed, let A be a sequence for which this criterion holds and assume 
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that we already integrated out the first i — 1 parameters. After mapping the integration domain of 
dx aj to [0, 1] we may assume that the integrand is of the form 

h = ££/;[wy,i|...|w ;> ], (4.2) 

m j 

where fj are algebraic functions and [to,- 1 \~-\CQj m ] G 33 m {Q. n ) for some m and n. The denominators 
of fj, CO/, i, 03j m are irreducible polynomials which map to the members of S/. If these polyno- 
mials are linear in x Gj then there is a k such that we can map these polynomials to the denominators 



of the 1 -forms in Q.^ of eq. 3.2 and we can express /; by terms of the form of eq. 3.15. 

It is an important advantage of this approach that by the polynomial reduction algorithm, i.e. 
by simple operations on polynomials and without integrating, we can decide whether the method 
applies and which order of parameters we should choose. Another advantage of the use 
is that if the polynomial reduction can be done, ^-linear combinations of these functions are 
sufficient to express the intermediate results after each integration. 

It was shown in reference [0] that the method is applicable for a large class of graphs and 



it is well-known how to relate vacuum-graphs to contributions to two-point functions []15[]. For 
graphs with further legs and with non-zero particle masses, we have to take the second Symanzik 
polynomial into account. Certain properties of this polynomial ^ give rise to the hope that the 
polynomial reduction and the above method can be extended to a large number of such Feynman 
graphs as well. It is furthermore important for us to move beyond the restriction of primitive graphs. 



Recent work of Kreimer and one of us [11, 12] provides a strategy to express Feynman integrals 



with UV-subdivergences by integrals for which the above method can be applied. The methods of 



[24] suggest that this method also generalizes to gauge theories. 
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